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1, INTRODUCTION 

Fault type and faulted phase identification are necessary in power distribution systems to keep 
uninterruptable power flow. Correct faulty phase identification is important to prevent switching of the 
incorrect phase or unnecessary three phase switching. It 1s also important in power system distribution 
network to quickly restore the system, improve system reliability, quickly repair the faulted line and reduce 
operating cost. Classifying the type of fault will help the relay to support several method components related 
with various fault events. Detecting the correct faulty phase also supports to satisfy single-pole tripping and 
auto reclosing requirements and to cut down the faulty phase(s). It may lead to increase the stability margin 
of the power system and possibly prevent needless loss of electricity in some regions. 

Most of fault classification methods mainly use various features of wavelet transform (WT) of 
voltage and current signals from the system measuring terminal. The fault type identification method in [1] 
classifies different fault types based on wavelet singular entropy (WSE) compared with the predefined 
threshold values. Local energy (LE) of wavelet transform of voltage waveform is applied in [2]. Moreover, 
the approximation coefficient-based fault type classification algorithm in [3] 1s considered based on the 
approximation coefficients of current traveling waves with the quarter cycle window and the method classify 
based on wavelet energy entropy (WEE) and wavelet entropy weight (WEW) [4-7], and the recursive wavelet 
transform method [8]. However, the outcomes of the algorithms in [1-8] are based on threshold values. Some 
fault type classification methods used the artificial intelligent methods, such as support vector machines 
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(SVM) [9-11], artificial neural network (ANN) [12], ANN with the use of particle swarm optimization (PSO) 
[13] and feedforward neural network combined with S-transform [14]. The fault classification algorithms 
identify various fault types based on sample values of voltage and current signals compared with their 
predefined values in [15] and time frequency characteristics of fault waveforms in [16]. Moreover, to 
enhance performance of directional relaying, fuzzy inference system is used to detect fault type with the help 
of magnitude of fundamental current and voltage signals [17]. Moreover, combined fuzzy logic and DWT 
based on the high system reactance to resistance ratio [18] and fuzzy logic based on prefault current and the 
variation of current [19] are proposed as well. However, the above artificial intelligent (AI) methods using 
wavelet coefficients are needed to train with the pre-estimated fault conditions. Morphological fault detection 
algorithms were also introduced in [20, 21] using arrival time, polarities and amplitudes of phase voltage and 
current signals, and polarities of DC component-based method [22]. But these algorithms are complicated to 
proceed and using polarity to detect faulty phase need to get the signals from both sides of the faulted line. 
Therefore, these algorithms cannot be used to identify by recording only sigle terminal of loop networks. 

The next fault identification methods firstly transform phase to modal components of voltage and 
current traveling waves and then extract fault features of these modal components with the use of discrete 
wavelet transform (DWT) filter. Among the various types of modal transform method, Karenbauer transform 
is used in [23, 24] and Clarke transform is used in [25]. The fault classification methods in [23-25] can be 
applied to identify various fault types for specified single transmission test system and sampling frequency. 
However, it cannot be used to detect fault types for multiple transmission systems and distribution networks 
because it used predefined threshold values for all fault types. The wavelet transform modulus maxima 
(WTMM) values may change according to fault types, faulted lines, fault distances and sampling frequency. 
These types of algorithms may not be used to cover for all types of variation which are fault resistance, 
faulted line, fault inception time and fault location from the measured terminal. 

In previous work [26, 27], the amplitude maxima of detail coefficients of modal current components 
are used to classify fault type and faulty phase by comparing the classification parameters each other 
according to basic power system fault principles. However, this method may identify wrongly when the 
magnitude and frequency of the original signals are very small because the detail coefficients represent the 
high frequency content of the signals. 

In this paper, a sample reflected current signal-based fault type classification algorithm 1s introduced 
using Karenbauer modal transform and time delay values of phase and modal current reflected signals. This 
technique utilizes only three-phase current signal from the measured terminal of the network. The approach is 
divided in to two parts. In the first step, the algorithm will detect the grounded or ungrounded short circuit 
faults based on the absolute maxima (AM) value of zero modal component. In next step, the algorithm will 
identify fault type and faulted phase under the grounded or ungrounded situation according to the decision of 
the previous step by using time delay values of phase and modal reflected current signals. 

The paper is organized as follows. Modal transform is presented in section 2. Wavelet transform 1s 
illustrated in section 3. Fault type identification is discussed in section 4 and section 5 presents simulation 
study for faulted phase classification. Section 6 display results and discussions and the last section is 
conclusion. 


2. MODAL TRANSFORM 

In this study, to identify the fault type, the phase current signals are needed to transform to modal 
current components to avoid mutual effect among them. All the distribution lines are expected to be fully 
transposed, and the Karenbauer transformation matrix is applied as a modal transformation matrix which 
is follows: 


Tq 1-1 0 
Ig} _ 111 0 -1 ‘a 

—- v4 I l 
Iy| 3111 «21 a (1) 
ie Oa ale 


where a,b,c stand for phase A, B, C respectively and a,8,0,y stand for the modal components. The modal 
transformation coefficients are real numbers and can be inferred by three phase instantaneous sampling data. 
This characteristic is very useful for transient signal analysis. The matrix of Karenbauer transform is a full- 
order matrix. o-modal refers to the line modal between phase “a” and phase “b’, B-modal refers to the line 
modal between phase “a” and phase “c”’, O-modal refers to the residual current component of three-phase 
current signals and y-modal refers to the line modal between phase “b” and phase “c’’. 
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3. DISCRETE WAVELET TRANSFORM 

The discrete wavelet transforms (DWT) is an effective time-frequency signal-processing tool which 
can be used to analyze dynamic sampled signals. In past decades, the DWT has been commonly utilized to 
analyze the dynamic oscillation of power system disturbances. 

Basically, the DWT separates the frequency-band of the input signal into low- and high-frequency 
components, which are called here as approximation and detail coefficients, respectively. In this paper, 
only arriving times of detail coefficient are used to detect fault type and faulty phase. Given a digital signal 
function x(k), its coefficients can be computed as follows: 


DWT(m,n) = Teeny (Se) 2) 


where do is the scale factor, bo is the translation factor and k is integer variable and it refers to a sample 
number of an input signal. 


4. TRAVELING WAVE THEORY 

In the proposed method, the recorded three-phase current waves from the substation are transformed 
to modal components. The key of this paper is to utilize the arriving time delay between the phases and 
modal current reflected signals come from fault point. Time delay between zero mode current and phase 
currents formula in (3) and time delay between zero mode current and modal currents formula in (4) are 
given as follows: 





At, =to-t, p=a.b,c, (3) 
Ati, —to tn sm=a,B,y, (4) 
je = (5) 
Vo 
xX 
t= =, 6 
p ata => (6) 
1 
Vo), (7) 
LoCo 
1 
i ————— (8) 


VLC, 


where At, and At,, are the time delay values of phase and modal current reflected signals, to is the arriving 
time of zero mode current traveling wave, t, and t,, are the arriving times of phase and modal current 


traveling wave signals. The arriving times can be calculated from vy, v, and x shown in (5) and (6). V, and 
v, are the velocities of the aerial and zero mode components of traveling waves respectively, and xis the 
distance between fault position and measured terminal. These velocities are calculated from characteristic 
impedance of the distribution line shown in (7) and (8) respectively. Ly) and Cy, are zero sequence inductance 
and capacitance of the power line and L, and C, are positive sequence impedance of the power line. When a 


fault occurs in the system, the characteristic impedance will change according to fault section and faulted 
phase. At that time, the traveling wave velocities will vary related to characteriatic impedance and the 
arriving times of all components will differ according to faulty phase. Therefore, the proposed technique used 
varing patterns of time delay values come from different arriving times related to fault type and faulty phase. 
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5. FAULT TYPE IDENTIFICATION 

The proposed technique is divided in to two parts. In the first part, the algorithm will detect the 
grounded or ungrounded short circuit faults based on the absolute maxima (AM) value of detail coefficient of 
zero modal component. The algorithm will finally identify fault type and faulted phase under the grounded or 
ungrounded situation according to the decision of the previous step by using time delay values of DWT of 
phase and modal current reflected signals. 


5.1. Grounded or Ungrounded Fault Discrimination 

After recording the fault current signals, the three-phase currents are transformed to the modal 
components. Next, AM value of detail coefficient of zero model components is applied to decide whether the 
system is grounded or ungrounded fault. Figure | illustrates the flow chart for detecting the grounded and 
ungrounded faults by zero sequence component. 

For appropriate discrimination between grounded and ungrounded conditions, a threshold value e 
is proposed. This threshold value represents the absolute maxima value of detail coefficient of zero modal 
because according to theory, the residual current can have under system unbalanced conditions, on the other 
hand, this value may be zero under system balanced condition. The threshold value should be assumed based 
on applied system situation and sampling frequency. Therefore, for the test system in this paper according to 


simulation results, it 1s found that for grounded cases the value is over 1x10° except for three phase-to- 
ground faults. In contrast, for ungrounded faults these values are below 1x10’. Therefore, it is chosen to 
define a definite threshold value (¢) as 1x10°' to cover numerical errors for accurate decision between the two 
types. In the following flow chart, 





dy| is absolute maxima (AM) value of detail coefficient of zero 


sequence current. 
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Figure 1. Fault type identification flow chart 













































































5.2. Classification of Fault Type 

The faulty phase is identified by utilizing six parameters. They are time delay values of phase and 
modal current reflected signals, Dt, ,Dt,,Dt,,Dt,,Dt,andDt, receptively. Fault type can be determined 
using the following procedure: 
a. Record current signals at main substation (1,, Ip, I.) 


b. Transform phase current components to modal components 
c. Apply Discrete Wavelet transform (DWT) with db4, db6 and db8 at level | and 
d. Calculate time delay values (Dt, ,Dt,,Dt,,Dt,,Dt,andDt, ) of detail coefficients of all phase and 
modal current signals. 
For grounded fault, firstly the algorithm will check whether system fault is single line to ground 
(SLG) or two-phase-to-ground (LLG) by using the time delay values of modal current reflected signals 
(Dt, ,Dt,andDt, ). The time delay is the time difference between the arriving times of zero mode current 
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signal and phase or modal current reflected signals from fault point to the measured substation. If system 
fault is LLG fault, the time delay values of modal current components must be the same because every modal 
component are related with one or two faulted phases. In that case, the faulted-phase identification algorithm 
will detect faulted phases by using the time delay values of phase reflected signals (Dt,,Dt,andDt. ) and the 


arriving time of faulted-phase signals may differ from that of the healthy signal. Otherwise, if system fault is 
SLG fault, the time delay values of a modal component related with two healthy phases may differ from that 
of other two modal components which are related with this faulted phase. For ungrounded fault, the proposed 
ungrounded faulted-phase identification algorithm will detect faulted-phase by using the same procedure 
used in LLG fault. 

The grounded faults identification flow charts are described in Figure 2, 3 and 4. In the first step, 
the algorithm will check whether the fault is SLG or LLG by comparing sum of absolute maxima values of 
modal components and residual current because the high frequency content of residual current must less than 
that of sum of all modal components when LLG fault occurs in the system, and this identification procedure 
is described in Figure 2. 

When SLG fault occurs in the system, the classification algorithm displayed in Figure 3 will activate 
and continuously select the faulted phase. In this process, the modal component related only healthy phases 
have minimum time delay values because of variation of traveling wave velocity. When a disturbance 
happens in the system, the traveling wave velocity will vary according to system characteristic impedance via 
fault assisted. For example, when a SLG (phase A to G) fault occurs in the system, the arriving time of the 
reflected signal of the faulted phase will arrive later than that of healthy phases depend on fault resistance. 
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Figure 2. Grounded fault type identification Figure 3. Single line to ground faulted phase 
flow chart identification flow chart 


Similarly, the two-phase to ground fault procedure is illustrated in Figure 4. The arriving time of 
reflected signal of healthy phase may arrive earlier than that of faulty phases. For ungrounded fault, 
the algorithm will discriminate directly, and the process 1s shown in Figure 5. If phase to phase fault occurs, 
the time delay value of the healthy phase current may have less than that of faulted phases because the 
arriving times of faulted phase will later than that of healthy phase because the traveling wave velocities of 
healthy and faulted-phase are different as system characteristic impedance is different due to fault. The time 
delay can be calculated from the arriving times of zero mode and phase current reflected signals. In this 
algorithm, 4 is the predefined values of time delay difference between modal current reflected signals. In this 
paper, the value of 1, 1x10* is assumed because the algorithm will use 4 to determine whether the time 
delay difference between modal current signals is zero or not. Therefore, if the time delay diference between 
modal current signals is zero, the time delay difference values at the specified test cases must greater than 1. 
The simulation results and performance of all the proposed algorithms can be seen in the following sections. 
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Figure 4. Two line to ground faulted-phase Figure 5. Ungrounded fault type identification flow chart 


identification flow chart 


6. SIMULATION RESULTS AND DISCUSSIONS 

In this paper, evaluation of fault type identification is evaluated on a power system simulation as 
demonstrated in Figure 6. The all fault signals of the loop distribution system is generated by using 
MATLAB SIMULINK program. 


B-T 


B-2 





Figure 6. IEEE 14-Bus modified test system 


The test network is supported by three power supplies of 33kV/10kV, 50Hz and a loop overhead 
distribution system (IEEE 14-bus modified test system). The faults are tested on 16 points on all branches of 
the test system. All fault types and faulty phases are tested with various fault impedances (R¢) (0.001Q-50Q) 
and different inception time (0.06sec, 0.02sec, 0.0005sec) of faults are studied. For every event, the time 
delay values of phase and modal current reflected signals with mother wavelets (db4, db6 and db8) at level 1 
are applied to discriminate fault type and faulted phase. It is found that these values are different for all fault 
types (SLG, LL, LLG, LLL and LLLG) with different impedances, feeders, locations and inception times. 
However, maintaining their characteristics of basic theory can be found. The simulation time is 0.1s with 
time step 0.lus. System parameters are generator source resistance 13.96Q and source inductance of 0.35H. 
Distributed line parameter model is used as overhead distribution line parameters. 


Indonesian J Elec Eng & Comp Sci, Vol. 17, No. 2, February 2020: 551 - 563 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 557 


6.1. Determination of Fault Type 

In this section, the fault type and faulty phase are discriminated by using time delay values of phase 
and modal current reflected signals. The estimated time dealy values are as shown in Figure 7 and 8 as an 
example. In this example, when phase ‘c’ to ground fault happens tn the test system, the amplitude maxima 
of zero sequence current is 0.001829 greater than threshold value (¢€). According to the algorithm shown in 
Figure 1, system fault is grounded fault. In that case, the time delay values are Dt, =9.73' 10 = 
Dt, =Dt, = 0.0001242, and these time delay values of modal current signls are not the same and the time 


delay of phase current are the same. Therefore, fault is SLG fault. For selecting correct faulted phase, 
the proposed algorithm will determine that the phase with the minium or maximum time delay is as the 
faulted phase. For this example, system fault is phase ‘c’ to ground fault because Dt, is different from 


Dt, andDt, . Therefore, fault type is single line to ground fault (phase ‘c’ to ground). 


The second example is phase ‘b’ to phase ‘c’ (BC) fault shown in Figure 9. When phase ‘b’ to phase 
‘c’ fault occurs in the test system, the phase current time delay values are At, = 0.0003129, 


At, =At, =0.0003665 because the time delay values of modal current signals must be the same and that of 


faulted-phase current signals must be different from healthy phase current signals according to traveling 
wave theory. Therefore, the arriving time of phase ‘a’ differs from that of other two phases because the 
traveling wave velocities of healthy and faulted phase are different as system characteristic impedance is 
different due to fault. This condition means system fault 1s phase “b’to phase ‘c’ fault. Other fault types can 
be identified as this example and test performance of proposed algorithms for various fault types are shown 
in the following sections. 
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Figure 7. The arriving time of modal current signals under phase ‘c’ to ground fault condition 
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Figure 8. The time delay value of reflected modal current signals under phase ‘c’ to ground fault condition by 
different mother wavelets (1) db4 and (11) db6 
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Figure 9. The time delay value of reflected phase current signals under BC fault condition by different 
mother wavelets (1) db4 and (11) db6 


6.2. Testing Results of the Proposed Algorithms 

The proposed method is tested with various parameters and shows high accuracy and meeting the 
condition set for each case with different faulty feeder, different fault resistance (R¢), different fault location 
along the same faulted line and inception time, as illustrated in Tables | to 4. 

In Table 1, the test results of proposed fault identification algorithm can be seen with minimum and 
maximum fault resistance conditions. According to this table, the content of high frequency transient can be 
varied by fault resistance. However, the proposed method can identify correctly because it detects the fault 
type and faulted phase by comparing the detection parameters each other according to basic principle of 
power system fault. 

In Table 2, the test results of grounded and ungrounded faults are illustrated according to different 
faulted feeders. The faulted feeder can affect to vary the time delay values of fault type identification 
parameters because the three-phase current signals are measured from only one substation of the loop 
distribution grid and the arriving times of faulted and healthy phases can be varied depend on the faulted line 
of the distribution network. 

Similarly, in Table 3 and 4, performance of proposed method with various fault inception times and 
different fault locations are illustrated respectively. In Table 3, the variation of fault type detection 
parameters also depends on the fault inception time because if a fault occurs at a time while the original 
signal crossing zero point (40 milli-sec and 60 milli-sec), there have only the effect of fault transient. 
Otherwise, when a fault occurs at other non-zero point (2.5 milli-sec), the delay time of reflected signal can 
different depend on the fault inception time. Moreover, according to Table 4, the time delay of phase and 
modal reflected signals are different according to the fault distance away from the measured terminal. Even 
though fault occur at the same faulted line and distance, the time delay values of faulted and healthy phases 
are not the same because their characteristic impedances are the same when a fault happens in the system. In 
this study, fault occur at the line L7, the fault distance is 1 km, 2 km and 3 km away from bus B-5, the lowest 
time delay values can be got when fault occur at the closest distance (1 km) from the measured terminal. 

The proposed algorithm discriminate fault type and faulty phase by comparing the time delay values 
of all modal and phase reflected signals each other without using a predefined value. Therefore, the proposed 
method can be used to classify the correct fault type and select faulty phase under various situations, such as 
various fault resistances, all faulted lines, various inception times and different fault locations along the same 
faulted line. 


Table 1. Performance of Proposed Fault Type Identification Algorithm for all Fault Types with Rf=0.001 
ohms and 50 ohms 


Fault Re Diy Dip Dig Dtg Dtp Di, 

T Oh Result 
ype (Ohm) 103 ‘10 3 47 3 40 3 rao ‘40 3 

CG 0.001 0.0971 0.1241 0.1241 1.181 1.181 1.181 CG 
CG 50 0.0971 0.1241 0.1241 1.181 1.181 1.181 CG 


BCG 0.001 0.1241 0.1241 0.1241 0.3129 0.3665 0.3665 BCG 
BCG 50 0.1241 0.1241 0.1241 0.3129 0.3665 0.3665 BCG 
BC 0.001 0.3665 0.3665 0.3665 0.3129 0.3665 0.3665 BC 
BC 50 0.3665 0.3665 0.3665 0.3129 0.3665 (0.3665 BC 
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Table 2. Performance of Proposed Fault Type Identification Algorithm for SLG (AG) 
Fault with R=0.001 ohms for All Branches 


Fault Faulted 
Type Line 
AG Ll 
AG L2 
AG L3 
AG L4 
AG | Bes) 
AG L6 
AG LZ 
AG L8 
AG L9 
AG L10 
AG L11 
AG | a bes 
AG L13 
AG L14 
AG L15 
AG L16 


Dt 


a 

“10 3 
0.000 
0.000 
0.0989 
0.7030 
0.0586 
0.0937 
0.1241 
0.0821 
0.1173 
0.1173 
0.1173 
0.2245 
0.2245 
0.1173 
0.1591 
0.2403 


Dtp 


"10 3 
0.000 
0.000 
0.0989 
0.7030 
0.0586 
0.0937 
0.1241 
0.0821 
0.1173 
0.1173 
0.1173 
0.2245 
0.2245 
0.1173 
0.1591 
0.2403 


Dt, 


10 3 
3.300 
3.00 
0.1054 
0.3040 
0.0187 
0.0537 
0.0971 
0.0555 
0.0839 
0.0905 
0.0905 
0.1711 
0.1711 
0.0905 
0.1323 
0.2005 


Dt 


a 

10 3 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 


Dtp 


10 3 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 


Dt 


Cc 
"10 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 


Result 


AG 
AG 
AG 
AG 
AG 


Table 3. Performance of Proposed Fault Type Identification Algorithm for All Fault Types with Rs=50 ohms 
for Various Fault Inception Time 


Fault pant 
Type Inception 
Time (sec) 
AG 0.06 
AG 0.04 
AG 0.0205 
BCG 0.06 
BCG 0.04 
BCG 0.0205 
BC 0.06 
BC 0.04 
BC 0.0205 


lq 

(ms) 
0.1241 
0.1241 
0.1241 
0.1241 
0.1241 
0.1241 
0.0939 
0.3665 
0.0939 


Dtp 


(ms) 
0.1241 
0.1241 
0.1241 
0.1241 
0.1241 
0.1241 
0.0939 
0.3665 
0.0939 


Dt, 


(ms) 
0.0971 
0.0971 
0.0971 
0.1241 
0.1241 
0.1241 
0.0939 
0.3665 
0.0939 


Dt, 

(ms) 

1.181 
0.0001 
0.0001 
0.0971 
0.3129 
0.0001 
0.0069 
0.3129 
0.0469 


Dtp 

(ms) 

1.181 
0.0001 
0.0001 
0.1241 
0.3665 
0.1241 
0.0939 
0.3665 
0.0939 


Dt 

(ms) 

1.181 
0.0001 
0.0001 
0.1241 
0.3665 
0.1241 
0.0939 
0.3665 
0.0939 


Result 


Table 4. Test Results of Proposed Fault Type Identification Algorithm for All Fault Types with Re=50 ohms 


for Different Fault Distances Along the Line (L7) 


Fault Fault 
Type Distance 
AG lkm 
AG 2km 
AG 3km 
BCG lkm 
BCG 2km 
BCG 3km 
BC lkm 
BC 2km 
BC 3km 
ABC 1km 
ABC 2km 
ABC 3km 


6.3. Filter Analysis 


Dt, 


(ms) 
0.0821 
0.1241 
0.1056 
0.0821 
0.1241 
0.0787 
0.0819 
0.3665 
0.1055 
0.0555 
0.0973 
0.0555 


Dtp 


(ms) 
0.0821 
0.1241 
0.0787 
0.0821 
0.1241 
0.0787 
0.0819 
0.3665 
0.1055 
0.0555 
0.0973 
0.0555 


Dt, 


(ms) 
0.0555 
0.0971 
0.0787 
0.0821 
0.1241 
0.1056 
0.0819 
0.3665 
0.1055 
0.0819 
0.1173 
0.0819 


Dt, 


(ms) 
0.0001 
0.0001 
0.0001 
0.0555 
0.3129 
0.0787 
0.0555 
0.3129 
0.0785 
0.0555 
0.0973 
0.0555 


Dtp 


(ms) 
0.0001 
0.0001 
0.0001 
0.0821 
0.3665 
0.1056 
0.0819 
0.3665 
0.1055 
0.0819 
0.1173 
0.0819 


Dt. 


(ms) 
0.0001 
0.0001 
0.0001 
0.0821 
0.3665 
0.0787 
0.0819 
0.3665 
0.1055 
0.0819 
0.1173 
0.0819 


Result 


AG 
AG 
AG 
BCG 
BCG 
BCG 
BC 
BC 
BC 
ABC 
ABC 
ABC 


In this section, level 1 detail coefficients of the three mother wavelets (db4, db6 and db&8) are 
applied to detect correct fault type and faulty phase, and then compare their outputs and choose the suitable 
filter with the highest accuracy because the arriving times of reflected signals passed through these three 
filters may different as these three mother wavelets have different number of non-zero coefficients and 
different amplitude of non-zero corfficients. The following tables are the comparison outcome of the three 
mother wavelets with respect to different conditions, such as various fault resistance, different fault inception 
times, different faulted lines and different fault distance from the measured terminal. 

In Table 5, the level 1 detail coefficients of the phase and modal current reflected signals by three 
mother wavelets are used to classify fault type and faulted-phase for all lines (L1 to L16) of test system by 
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creating single line to ground (SLG), phase ‘a’ to ground (AG) faults. According to the result of Table 5, 
the accuracy of db6 1s the highest accuracy in those three mother wavelets (db4, db6 and db8). 

In Table 6, there are eight cases for all fault types with respect to fault resisatnces (0.001 Wand 
5O0W) by using the detail coefficients of three filters (db4, db6 and db&8) respectively. In that case, 
the proposed algorithm can classify all fault types accuratly and the average fault type and faulted phase 
identification accuracy for all filters are 100 percent because the value of fault resistance cannot influence on 
the arriving time via the time delay of phase and modal current reflected signals. 


Table 5. Performance of Proposed Fault Type Identification Algorithm for SLG (AG) Fault by using db4, 
db6 and db8 with R=0.001 ohms for Different Faulted Lines 


Case Actual Faulted Estimated Fault Type by 
Fault Type Line db4 db6 
1 AG L; AG AG AG 
2 AG L» AG AG AG 
3 AG L3 AG AG CG 
4 AG Ly CG AG AG 
5 AG Ls AG AG CG 
6 AG Le AG AG CG 
i) AG Ly AG AG AG 
8 AG Ls AG AG AG 
9 AG Lo CG AG AG 
10 AG Lio CG AG AG 
11 AG Li, AG AG AG 
lip AG Liz AG AG AG 
13 AG Lis AG AG CG 
14 AG Li4 CG AG AG 
15 AG Lis AG AG AG 
16 AG Lie AG AG CG 
Average accuracy (%) 1 100 68.75 


Table 6. Performance of Proposed Fault Type Identification Algorithm for All Fault Types by using db4, 
db6 and db8 with Different Faulted Resistances at L7 


Actual Fault resistance Estimated Fault Type by 
ack Fault Type (W) db4 db6 db8 
1 AG 0.001 AG AG AG 
2 AG 50 AG AG AG 
3 AG 100 AG AG AG 
4 ABG 0.001 ABG ABG ABG 
5 ABG 50 ABG ABG ABG 
6 ABG 100 ABG ABG ABG 
2 AB 0.001 AB AB AB 
8 AB 50 AB AB AB 
9 AB 100 AB AB AB 
Average accuracy (%) 100 100 100 


The performance of fault type classification with respect to fault inception time (0.0205, 0.04 and 
0.06 sec) by using three filters 1s shown in Table 7. In this case, all the test cases can be detected by three 
mother wavelets. Similarly, in Table 8. illustrate that the method can select the fault type and faulty phase 
according to fault distances (1 km, 2 km and 3 km of the same faulted line). Although difference locations 
can get different arriving times, the patterns of time delays for three filters (db4, db6 and db8) are the same 
for this case. Therefore, the average accuracy 1s 100%. 

According to Table 5 to 8, the mother wavelet db6 is the most suitable filter for the proposed 
technique because it has the highest fault type identification accuracy in different faulted lines shown in 
Table 5, the equal accuracy in other conditions are shown in Table 6, 7and 8, respectively. 
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Table 7. Performance of Proposed Fault Type Identification Algorithm for All Fault Types by using db4, 
db6 and db8 with Different Fault Inception Times at L7 


Ca Actual Fault Inception Estimated Fault Type by 
se Fault Type Time (s) db4 db6 db8 
1 BG 0.0205 BG BG BG 
2 BG 0.04 BG BG BG 
3 BG 0.06 BG BG BG 
A BCG 0.0205 BCG BCG BCG 
5 BCG 0.02 BCG BCG BCG 
6 BCG 0.06 BCG BCG BCG 
7 BC 0.0205 BC BC BC 
8 BC 0.04 BC BC BC 
9 BC 0.06 BC BC BC 
Average accuracy (%) 100 100 100 


Table 8. Performance of Proposed Fault Type Identification Algorithm for All Fault Types by using db4, 
db6 and db8 with Different Fault Distances at L7 


Case Actual Fault Distance Estimated Fault Type by 
Fault Type (km) db4 db6 db8 
1 BG 1 BG BG BG 
2 BG 2 BG BG BG 
3 BG 3 BG BG BG 
4 ABG 1 ABG ABG ABG 
5 ABG 2 ABG ABG ABG 
6 ABG 3 ABG ABG ABG 
i AB 1 AB AB AB 
8 AB 2 AB AB AB 
9 AB 3 AB AB AB 
10 ABC 1 ABC ABC ABC 
1] ABC 2 ABC ABC ABC 
12 ABC 3 ABC ABC ABC 
Average accuracy (%) 100 100 100 


7. CONCLUSION 

This paper proposes a simple fault identification technique in loop distribution systems. 
The proposed technique passed through Karenbauer Transform and Discrete Wavelet Transform filter to 
estimate a vigorous result for determination of fault type and faulty phase by using only current traveling 
waves at one source bus of the test distribution network. In all algorithms, the time delay values of modal 


current reflected signals and phase current reflected signals (Dt, ,Dt, ,Dt 2 ,Dt,,Dt, and Dt.) are applied as 


the fault classification parameters to identify fault type and faulty phase in loop distribution systems which 
are very complex. 

The proposed method classifies the fault type and faulty phase correctly by comparing these 
classification parameters to each other instead of using threshold values except for identifying whether 
system fault is grounded or ungrounded case. If threshold values are used for the fault classification process, 
it need to be defined by evaluating for every test system. However, for the proposed method, the fault type 
and faulty phase can be classified correctly without defining and evaluating threshold values for every cases 
and the three mother wavelets (such as db4, db6 and db8 which are suitable for high frequency power system 
signals) are used to ensure and improve the reliability of the proposed method. The integrity of the proposed 
technique is successfully tested with data obtained by MATLAB simulations with various cases of fault 
current signals at different conditions (faulted feeders, faulted resistance, fault locations and inception times). 
This proposed technique can be used to implement in real data with actual fault records. 
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